Quantum communication cost of preparing multipartite entanglement 
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We study the preparation and distribution of high-fidelity multi-party entangled states via noisy 
channels and operations. In the particular case of GHZ and cluster states, we study different 
strategies using bipartite or multipartite purification protocols. The most efficient strategy depends 
on the target fidelity one wishes to achieve and on the quality of transmission channel and local 
operations. We show the existence of a crossing point beyond which the strategy making use of 
the purification of the state as a whole is more efficient than a strategy in which pairs are purificd 
before they are connected to the final state. We also study the efficiency of intermediate strategies, 
including sequences of purification and connection. We show that a multipartite strategy is to be 
used if one wishes to achieve high fidelity, whereas a bipartite strategy gives a better yield for low 
target fidelity. 

PACS numbers: 03.67.Mn, 03.67.Hk, 03.67.Pp 



I. INTRODUCTION 

In the past, (multipartite) entanglement has been 
mainly considered as a puzzling artifact of quantum me- 
chanics. More recently, however, the focus on entangle- 
ment has shifted, as it was realized that entanglement 
also constitutes a valuable resource for quantum infor- 
mation processing. Possible applications of multipartite 
entanglement include certain security tasks in distributed 
communication scenarios 0, Q , the improvement of fre- 
quency standards Q , as well as measurement based quan- 
tum computation schemes Q. 

In this context, the problem of generating multipar- 
tite entanglement of high fidelity arises. If entangled 
states are to be distributed among spatially separated 
parties, as it is e. g. required in distributed communi- 
cation scenarios, the main obstacle comes from chan- 
nel noise. Possible ways to overcome channel noise and 
hence to successfully generate high-fidelity multipartite 
entangled states have been developed. These methods 
are based on (i) quantum error correction and make use 
of (concatenated) quantum error correction codes [jj, or 
(ii) entanglement purification [S El El ■ While (i) is 
applicable to directly distribute arbitrary states, (ii) con- 
centrates on the generation of specific, maximally entan- 
gled pure states. The generation of maximally entangled 
pairs of partides allows in turn to distribute arbitrary 
states by means of teleportation. In both cases, a sub- 
stantial overhead is required to guarantee successful, high 
fidelity generation of the desired states. In (i) this over- 
head arises from redundant encoding, enabling one to 
perform error correction, while for (ii) several identical 
copies need to be prepared and locally processed to gen- 
erate high fidelity entangled states. The quantification 
of this overhead, or the quantum communication cost, 
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which we shall define more precisely, is the main concern 
of this article. 

To be specific, we will concentrate on schemes based on 
entanglement purification. Thcse schemes are specially 
suited to generate entangled states of a specific form, 
and are hence expected to perform better than general 
purpose schemes such as (i). In fact, a remarkable ro- 
bustness of entanglement purification protocols against 
noise in local operations -which we consider in addition 
to channel noise- has been found [ÜE!- That is, errors of 
the order of several percent in local control operations can 
be tolerated, still allowing for the generation of high fi- 
delity entangled states, even in the prcscncc of very noisy 
quantum channels and with only a moderate overhead. 
For perfect local operations, the required overhead in re- 
sources is solely determined by the noisy quantum chan- 
nels. In this case, the channel capacity [lfj, [ll| provides 
a suitable measure for this overhead. In a bipartite com- 
munication scenario, the channel capacity gives the op- 
timal rate of quantum communication, i. e. the amount 
of quantum information transmitted per actual channel 
usage. While one might think that the abstract notion of 
channel capacity may also be employed to our problem 
-the generation of certain high fidelity entangled pure 
states-, one actually faces a number of difficulties. First, 
channel capacities are asymptotic quantities which are 
very complicated to calculate; second, the definition of 
channel capacity is not suitable to account for imperfect 
local operations (e. g. noise in local coding and decod- 
ing procedures); and third, we are actually considering 
a restricted problem, namely the generation of specific 
multipartite entangled states, rather than the successful 
transmission of arbitrary quantum information. 

We thus introduce a quantity closely related to quan- 
tum channel capacity, namely the quantum communica- 
tion cost Cf,g- Cf,g denotes a family of quantities which 
specify the number of uses of the noisy quantum chan- 
nel required to prepare a specific (multipartite) entan- 
gled state \G) with fidelity F > F. In this paper, we 
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will focus on target states \G) which are so-called two- 
colorable graph states. These states include, for instance, 
GHZ states and cluster states -a universal resource for 
measurement based quantum computation Q- and they 
are locally equivalent to codewords of Calderbank-Shor- 
Steane error correcting codes 0, El . We establish up- 
per bounds on Cf,g by optimizing over a large class of 
different strategies that generate these multipartite en- 
tangled states. These strategies include, as extremal 
cases, (i) the generation and purification of pairwise en- 
tanglcment, from which, by suitable conncction processes 
(or, alternatively, teleportation) the desired multipartite 
states are generated; (ii) the generation and direct mul- 
tipartite purification of the desired target states. Inter- 
mediate strategies, e. g. the purification of smaller states 
to high fidelity and their subsequent connection to the 
desired larger state, will also be investigated. Depend- 
ing on the actual noise parameters for channels and local 
control operations, and on the desired target fidelity F, 
the optimal strategy varies. For high target fidelities, 
multipartite strategies turn out to be favorable. 

This article is organized as follows: In Sec. |üJ we 
present the concepts we will use: We start with a review 
of the graph state formalism in ordcr to introduce nota- 
tion and the two types of states we wish to distribute, 
namely GHZ states and 1D cluster states. We shall also 
introduce a technique to connect two smaller graph states 
to obtain a larger one. Then, we give details for our noise 
models and review the employed purification protocols. 
Readers familiar with these concepts may skip this sec- 
tion. Sec. II I II explains the different strategies for em- 
ploying the protocols that we wish to compare. The ac- 
tual comparison is done using extensive numerical Montc 
Cario simulations and results are presented in Sec. II VI In 
order to corroborate these results we have done analyt- 
ical studies for certain restricted noise models ÍSec. Ivll. 
We conclude with a summary ÍSec. IVIfl . 



II. BASIC CONCEPTS 
A. The graph-states formalism 

A graph G — (V, E) is a collection V — {a, b, c, . . . } 
of N = IT^I vertices connected by edges E C [V] 2 . The 
description of the edges is given by the adjaceney matrix 
Te associated with the graph 

{1, if a and b are connected by an edge, 
i.e. {a,b}eE 
0, otherwise 

The neighborhood N a c V of vèrtex a is defined as the 
set of vertices connected with it by an edge, N a = {b : 
{a,b}£E}. 

With each graph G we associate a pure quantum state. 
If the graph's vèrtex set can be separated into two sets 
A and B such that no edges exist between vertices of 



the same set, we call it a two-colorable graph (Footnote 
[Ï4|). The vertices are qubits and the edges represent 
interactions. 

There are three equivalent descriptions of graph states 
which are reviewed in the following sections (For a de- 
tailed treatment see [l5jt : 



1. Graph states in the stabilizer formalism 
Associated with a graph G is a set of N operators 

4 a3 =4 o) n 

b£N a 

They form a complete set of commuting observables for 
the system of qubits associated with the graph and there- 
fore possess a set of common eigenstates which form a 
basis of the Hilbert space. These eigenstates are called 
graph states and are here written as \G,fi) where the 
a th component of vector fj, S {0, 1}^ is equal to if 

K^lCn) = \G,fi) and 1 ií K^lCfj.) = -\G,p). We 
abbreviate \G) :— \G, 0). We also sometimes suppress 
the letter "G" and write just if the context makes 
clear which graph G is meant. 



2. Graph states in the interaction picture 

A graph state with \x = can be written in the com- 
putational basis in the following manner: 

\G) = ( Yi Aziab) ) (2) 

\{a,6}e£ / 

where AZ is the controlled phase gate, 

AZ^ = |00) Qb (00| + |01) afc (01| + 

|10) o6 (10| - |11)„6<11|, (3) 

which corresponds to an Ising-type interaction, AZ^ ab ^ — 
e -iirH<· ab) ^ with interaction Hamiltonian i/( ab ) given by 

= |11)«6<H|. 

That is, \G) is generated from a pure product state by 

applying interactions between all pairs of partides con- 
nected by edges. 

We list some useful relations for later reference: The 
2 N common eigenstates of the operators can be 

generated from \G) by applying all possible produets of 
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FIG. 1: Producing a graph state in the VBS picture. 
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JV. This can be seen from 
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|G,0. 



which means that ai a) |G,0 
this relation, together with the fact that a 
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one can deduce the effect of ax an d (see Ref. |Ï5| 
and for proofs): Splitting the index vector /z into /z a , 
/xjv a (neighborhood of vèrtex a) , and fi Ra (remaining ver- 
tices), we write 

CJ i ^ ) \G,^aHN a HR a ) = \G,~p^HN a IJ,R a ) (5) 

a^\G,fj, a /j, Na n Ra ) = {-l) IJ ' a \G,fx a JïW2iJ,R a ) ( 6 ) 
o^ a) |G, n a fi Na fi Ra ) = i(-iy"-\G,-p^-fI]ï2n Ra ) (7) 
where the over-bar means bit complementation. 



3. Graph states in the valence bond sòlid (VBS) picture 

Another description of graph states was introduced in 
Ref. [lfjj. In this picture, every edge is replaced by a pair 
in a maximally entangled state, usually ( 1 00) + 1 01} + 
1 10) — |ll))/2. Each qubit a gets replaced by d a vir- 
tual qubits, where d a — \N a \ is the degree of vèrtex a. 
The physical qubit is recovered by projecting the virtual 
qubits onto the two-dimensional subspace of the physical 
one (see Fig. ^) using as projector 
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FIG. 2: Graphs for (a) cluster and (b) GHZ states. 



examples. By "cluster states", we mean graph states 
associated with a regular lattice as graph, in this article 
always a line as in Fig. |U(a), and with fi — 0. The term 
GHZ state will in this article be used for a graph state 
(again with fi — 0) associated with a star-shaped graph 
G* as in Fig. |2|(b). Such a state can be written as 



V2 



|0)® 



(JV-l) 



1 



\(g>(JV-l) 



and is hence in its entanglement properties equivalent 
to an "ordinary" GHZ state (10)®^ + (1)®^) = (1 <g> 

Had 8(JV_1) )|G») (where |±) = ^(|0) ± |1)) and Had is 
the Hadamard operation Had = -7= í J ^ j ) . 

5. Bell pairs and graph-state formalism 

In order to keep a certain homogeneity, we will employ 
a new notation for the states of the Bell basis, usually 
written as: 

I**} = ^=(|00)±|11)) 

I**) = ^ (|01>±|10>), 

Applying a Hadamard operation on the second qubit, 
one obtains a new basis formed by the graph states 
| G 2 , 00) , | G 2 , 01) , | G 2 , 10) and | G 2 , 1 1) , where G 2 denotes 
the graph composed of two vertices and one edge. Our 
new notation shows directly the relation between this two 
bases: 



|4>+) 
|*+) 
|$") 



\B;00) = Had (2) |G 2 ,00) 

\B;01) =Had (2) |G 2 ,01) 

\B;10) = Had (2) |G 2 ,10) 

IB; 11) = Had (2) |G 2 ,11). 



(9) 



P d = |Ò)<0...0| + |ï)(l...l|. 



(8) 



6. Connection of graph-states 



4- Cluster states and GHZ states 

In this article, we study the purification of graph states 
using two important representatives from this class as 



In this section, we define a procedure to connect two 
graph states, \G\) with N\ qubits, and |G 2 ) with N 2 
qubits, "fusing" together their respective vertices a\ and 
a 2 , yielding a state |G) with N\ + N 2 — 1 qubits. This 
process is depicted in Fig. [31 To realize this action, one 
applies a projective measurement on a± and a 2 , given by 
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n 2. Measure 

here in Y basis 
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3. Depending on measurement 
result, apply S or S f here. 



2 . n 



1 . Do local operations 
for neighborhood 
complementation 



-< 2. Connect 
1'» 3. Measure 1 /K S or S> 



4 \ \ 4. Correct for measure- 
2' • m' ment result 
3' 



P2 = |0>(00| + |1)(11| and = |0)(01| + |1)(10| (with 
outcomes and 1). P2 is defined like in the VBS picture. 
By similarity with this picture, if the result of the mea- 
surement is 0, the final state is the graph state resulting 
from the connection of G± and G2. If one obtains 1, a 
correction has to be done. As shown below, it is suffi- 
cient to apply IlbeN °* 1° ^ ne resulting state. Recall- 



ing that K { ç ] \G) = \G) with K t 

one sees that any graph state can be decomposed as 



(a) 



n 



b£N a 



\G) = |0) o ® |x> + |1> 
libeiv,, 



n 



ai^ to the state resulting from one obtains 



lx>- Applying 



FIG. 3: (Color online.) How to use the connection procedure 
of Sec. in A"6l to assemble large (a) cluster and (b) GHZ states 
from smaller cluster or GHZ states. 



n ai fc )p 2 ^ic7 1 )iG 2 )= n °í b) K d) io>®ixi>® n i»>+ii>® n \ = 



b£N ao 



b£N a 



deN a2 cGN ai 

|0) <8> Ixi) ® |X2> + 1] ^l 1 )® Ixi)® 1X2) - IG1+G2) (10) 

6e(JV ai +Ar a2 ) 



Fig. |3| shows how to use this technique to assemble 
cluster and GHZ states from IG2) states. 



B. Noise model 



Channel noise 



In any realistic setting, channels will be noisy. We 
study the influence of channel noise by considering re- 
stricted noise models, where the Kraus representation of 
the superoperators is diagonal in the Pauli basis. This 
is a common and usually sufficiently general model [Ït| 
(In particular, any noisy channel can be brought to such 
a form by means of (probabilistic) local operations). 
This allows for an efficient and convenient simulation by 
Monte Cario techniques (see Sec. II V À")l . We consider the 
following channels: 



depolarizing channel: 

p ^ £^{p) = qp + (4 a) p4 a) + 

a^pa^ + a^pa^ 
In case of depolarizing channel, we define 

p=(4«-l)/3 
which allow us to rewrite Eq. I|13|) as 
1-p 



£ {a) (p) 



pp 



4°V4 0) + °í a) p°í a) 



(13) 



(14) 



(1 — q) will be called the alteration probability and p the 
reliability. 



2. Local noise 



phase-flip channel: 

p^£^(p)=qp+(l-q)a^po 



(<>) 



(11) 



bit-flip channel: 

p^£Í a \p)=qp+{l-q)a^p<T^ (12) 



As part of the purification protocols, local one- and 
two-qubit unitary operations are employed which may be 
noisy. An imperfect operation is modeled by preceding 
the perfect operation í/( ab ) with the application of one of 
the noise superoperators £ from Eqs. I|11H13|1 . i. e. the 
state is transformed as 

p^U {ab) {£ a {£\p))) [/ t(Qb) . 
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We assume that the protocols are executed with the least 
possible number of operations to keep accumulatcd noisc 

low. Hence, if a two-qubit gate is preceded by one- 

qubit gates u[ a ^ and we apply one combined unitary 
jj(ab) = u[ a )u^u[f) w hich is subjected to noise only 
once. 



3. Commutation between connection and noise 

We now state an observation that will later (in Sec. 
IVBl be of use. 

For any graph states |Gi) and |G 2 ) which are con- 
nected by the procedure described in Sec. III A~6l one can 
show that the noise processes commute with the connec- 
tion procedure, if they are expressed by a superoperator 
by only a z Pauli operators. This comes from the fact 
that the neighborhood of the connected vertices 01 and 
a 2 changes with the connection and hence, a x and a y 
Pauli operators will affect different vertices (see Eqs. 101, 
© and 0). 

The commutation rules between projector P 2 (see Eq. 
iJHJ) and a z can be deduced from the following expression 
of the connected graph state: 

P 2 |G 1 )|G 2 )=P 2 (|0) ai |0) a2 | X }i|x)2 + 

+ LI ^í c) |0) ai |l) a2 |x)i|x)2 + 

c(EN a2 

+ LI ^ b) | 1 >a 1 |0) a2 |x)l|X>2 + 
b£N ai 

+ n ^ n ^ c) i i >a 1 ii> a2 ix>iix>2) 

beN ai c£N a2 

Recalling that <7 Z |0) = |0), <r z |l) = — 11) one can show 
that: 

P 2 a^\G x )\G 2 ) = a^P 2 \G x )\G 2 ) (15) 
P 2 a^\G x )\G 2 ) = a^P 2 \G x )\G 2 ) (16) 

C. Local noise equivalent 

To judge how close state p it to the desired state \tfj), 
one often usually employs the fidelity F := {ij)\p\tj)) . How- 
ever, is may be advantageous to regauge the fidelity by 
introducing the following derived measure: We define the 
local noise equivalent (LNE) as the level of local depolar- 
izing noise (in terms of the alteration probability (1 — q) 
of Eq. dU) that one has to apply to each qubit of the 
perfect state \tp) to deteriorate it to the same fidelity 
F as p has. The advantage of this measure is twofold: 
(i) It is more natural for uses of states in quantum er- 
ror correction schemes, as it can be compared directly 
to the fault-tolerance threshold in case of uncorrelated- 
noise models, (ii) It does not fall off exponentially with 



the size of a state for constant noise levels, as the fidelity 
does. On the other hand, it often cannot be calculated 
analytically in a straight-forward way We hence used a 
numerical Monte Cario simulation of the state deteriora- 
tion (which is why the LNE scale in the figures has error 
bars). 

D. Purification protocols 

The purpose of entanglement purification is the fol- 
lowing: One is given an ensemble of multi-party states, 
which all are distributed over two (or more) sites and ex- 
hibit entanglement between the sites. These states are 
only an approximation to the desired state | \1/ ) (\T/ 1 with 
an insufficient fidelity, which one wishes to improve. As 
the sites are spatially separated, one cannot apply joint 
operation on the distributed parts of a state. Instead, 
one compares (in case of the so-called recurrence proto- 
cols, which are considered here solely) pairs of entangled 
states, makes joint operations on them, and then mea- 
sures one of the state in order to gain information about 
the other. Only for specific measurement outeomes, the 
other state is kept. After iterating this procedure, one is 
left with an ensemble of smaller number of partides but 
highcr fidelity. 

1. Bipartite purification 

Several protocols have been proposed to purify bi- 
partite entangled states 0, 0, 0- To test the differ- 
ent strategies, we used the most efficient which can be 
used to purify an ensemble of |$+) states, namely the 
one described in Ref. 7|. We present here a modi- 
fied version of this bipartite entanglement purification 
protocol (BEPP) which allows for the purification of 
the connected graph-state pair. As we are concerned 
only with this graph in this section, we simply write 
\[iv) for the different basis states \G 2 ,\iv). Recali that 
|00) = 1/V2 (|Q)|+) + |l)|->) (see Eq. ©). 

Alice and Bob want to share entangled pairs with high 
fidelity. At the beginning they are given an ensemble of 
noisy 1 0, 0) states, each of them owning one part of the 
pairs. We consider a state diagonal in the graph-state 
basis, 

p = x o|00)(00|+x 01 |01)(01| + 

xi |10)(10|+x 11 |ll)(ll|. (17) 

We remark that such a Standard form can always be 
achieved by means of depolarization, i.e. applying cer- 
tain (random) local unitary operations. Each step of the 
protocol consists of the following operations: (i) Alice 
and Bob perform unitary operations on their partides, 
with Alice's and Bob's unitàries given by 
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(ii) Alice performs a CNOT operation from the first state 
to the second and Bob from the second state to the first; 

(iii) Alice and Bob measure the second state in different 
bases. To see the effect of this procedure, we calculate 
the fidelity and yield obtained after one step with two 
initial states given by ifTTj) . 

In (i), Alice and Bob apply S A and Sb, respectively, 
in order to swap |11) and |10). Then, in step (ii), they 
apply the bilateral CNOT. One can check that the effect 
of this operation on the graph state basis is given by the 
following map: 

\haHb)Wavb) >-> \ha © v a ,h b )\va,vb © Ms), (18) 

(Here, indicates bitwise and, i. e. addition modulo 2.) 
Last (iii), Alice and Bob measure the qubits of the target 
state. This is done in the eigenbasis {10)^, 1 1)^} of a x 
for Alice and in the computational basis {|0) z , |l) z }, for 
Bob. By this they obtain the eigenvalue of the correlation 
operator K 2 defined in Eq. (JIJ and determine the value 
of the second bit describing the state. If it is 0, they keep 
the first state. They discard it otherwise. 

After the measurement, they keep the control state 
with success probability k = (xoo + a^ii) 2 + (.xoi + xw) 2 
and the new coefficients are given by: 

• T oo = (^oo + x u)/k 

x oi = (^oi + a; io)/^ 
x' w = (2xooxu)/k 

x' u = (2x ix w )/k. (19) 

Hence, the fidelity is F = x' m = (xq + x\{)/k. The yield 
of the step, defined as the number of remaining states 
divided by the number of states before the step, is given 
by k/2 as half of the states (the targets) are measured 
and discarded. 

The unitary operations performed at the beginning of 
the protocol (step (i)) are required for its convergence. 
It guarantees that fidelity 1 is a fix point of the proto- 
col which is approached when iterating the procedure. 
The CNOT operation is a means of transferring infor- 
mation from the first qubit to the second. The mea- 
surement allows to distinguish between {|0, 0), |1, 0)} and 
{|0, 1), |1, 1)} and hence, determines the second bit of the 
index vector. 



2. Multipartite purification 

Multipartite purification protocols (MEPP) have been 
introduced in Re f. I l9l| for GHZ states, were further de- 
veloped in Ref. (23 and extended to all two-colorable 
graph states in Refs. [Hlï|. Recali that a two-colorable 
graph-state is a graph state in which the vertices can be 
separated into two sets Va and Vb such that no edges 
exist between vertices of the same set. Using the proce- 
dure described in Ref. Q, one can depolarize any mixed 



state p to the form 

P= V4,mJ G >/^,Mb)(G,/xa,Mb| (20) 

without changing the diagonal coefficients (where fiA, 
Hb are binary vectors corresponding to sets V A , Vb re- 
spectively). Hence we will restrict our attention to input 
states of this form. The protocol is composed of two 
subprotocols PI and P2 which we will describe here: 

a. Subprotocol PI: The states composing the en- 
semble of two-colorable graph-states are processed pair- 
wise. All parties belonging to set Va perform a CNOT 
operation from the second state of a pair of states to the 
first one while the parties belonging to set Vb perform a 
CNOT from the first one to the second one. This leads 
to the transformation 

\G,ha,Hb)\G,v a ,v b ) 

i-> \G, íia.íib ®vb)\G, va © ha,vb), (21) 

As in the bipartite protocol, the last step consists of mea- 
suring the second state of the pair. The parties belong- 
ing to set Va measure their qubit a in the eigenbasis 
{|0)xj |1)íb} 01 a x, obtaining results £ a S {0, 1}, while the 
ones belonging to set Vb make their measurement in the 
computational basis, obtaining results £{, £ {0, 1}. From 
this, we can calculate the part of the index vector of the 
measured state (second state of the r. h. s. of Eq. ijüj) 
corresponding to set V A : 

v' A = V A S VA = Ua © Cb ) 

V beN a I aeVa 

If this is 0, it is most likely that /j, a = 0, and hence, 
the first state is kept (and otherwise discarded). As con- 
sequence, in the expansion H2()(l of the ensemble density 
matrix, elements of the form Ao.^ B are increased. One 
finds that the new matrix elements are given by 

^~1A-~1B ~ ~ ^1A,VB^tA,y.B (22) 

where k is a normalization constant such that tr( j 5) = 1. 

b. Subprotocol P2: As explained above, subproto- 
col PI is employed to purify with respect to the eigenval- 
ues ha associated with set V A . The second subprotocol 
leads to the purification with respect to the eigenvalues 
of set Vb- It is obtained from PI by exchanging the roles 
of set Va and Vb ■ The protocol's action is described by 
the following map: 

\G, ha,Hb)\G, v>a,vb) 

^\G,HA®v A ,fi B )\G,u A ,VB®i'B)· (23) 

One measures the second state. The measurements on 
set Vb are done in the eigenbasis {10)^, |l) x } of a x while 
they are done in the computational basis in set V A . This 
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leads to the determination of part [13 of the index vector. 
As in subprotocol PI, one keeps the state ií Hb =0. The 
new coefhcient are given by 



(a) BPP strategy 



(b) MPP strategy 



{ (y A A ) | u a (Sha =1a } 



where k is a normalization constant such that tr(p) = 1, 
as before. 



III. STRATEGIES 



prepare and 
send pairs 



>M-y b 

purify pairs 



connect to obtain 
multipartite state 



locally prepare 
multipartite state 



send multipartite 
state through channels 



purify multipartite 
states 



A. Quantum communication cost Cf,g 

We now define our figure of mèrit, the quantum com- 
munication cost. We consider N spatially separated par- 
ties Ak, k = 1,2, ...,N which are pairwise connected 
by noisy quantum channels £ki, described by completely 
positive maps acting on density operators for qubits. 
We will quantify the quantum communication through 
these quantum channels using the quantum communica- 
tion cost Cki, i.e. the number of usages of the quantum 
channel £ k i , while classical communication between pairs 
of parties will be considered to be for free. Sending a sin- 
gle qubit through the quantum channel costs 1 unit, 
i.e. Cki = 1, while the transmission of an arbitrary state 
of M qubits costs Cui — M. We will be interested in the 
total quantum communication cost C, where 



(25) 



A'</ 



We consider the generation of multipartite entangled 
states (graph states, to be specific) \G) distributed among 
the parties A k . The goal is to generate states A = 
®\=\Pii where the fidelity of each p Í7 Fi = {G\pi\G), 
fulfills Fi > F. That is, each of the states has a fidelity 
larger than a threshold value F, which we call the "de- 
sired target fidelity". We remark that we demand that 
the ensemble of output states are in a tensor product 
form. In principle, weaker requirements such as that only 
the reduced density operators of A have fidelity larger 
than F are conceivable, however one faces certain diffi- 
culties in this case. For instance, it is not clear whether 
each of the copies of the state can be independently used 
for further quantum information processing tasks due to 
possible classical corrclations among the copies. Hence, 
we deliberately demand the tensor product structure. We 
will be interested in the total quantum communication 
cost C required to generate A = ®f =1 p% with Fi > F. 
In particular, we consider the quantum communication 
cost per copy, 



C 



F.G 



C 

T 



(26) 



where one optimizes over all possible strategies to gen- 
erate A. Due to this optimization, the quantity CV,g is 



FIG. 4: (Color online.) Distribution of iV-qubit GHZ states 
over noisy channels. (a) Bipartite entanglement purification 
strategy: Bell pairs are sent over the channels and purificd 
using a BEPP. The purified pairs are then connected (using 
the procedure of Fig. [3J to the desired GHZ state. (b) Mul- 
tipartite entanglement purification strategy: Alice prepares 
the GHZ state locally and sends all but one of the partides 
through the channels. Then, the MEPP protocol is used. 



very difficult to calculate. Hence we restrict ourselves to 
establish upper bounds on C'f,g by considering explicit 
strategies to generate high fidelity multipartite entangled 
states. 

Múltiple variations of this problem are conceivable. 
For simplicity we will assume that all parties are pairwise 
connected by identical quantum channels, £ = £ki- In- 
homogeneous situations where only some pairs of parties 
are connected by quantum channels (a restricted com- 
munication web), or where the classical communication 
is limited, or cases where quantum channels between dif- 
ferent pairs of parties are different (i.e. different noise 
parameter) will not be considered here. 

We will look mainly at two scenarios depicted in Fig. 21 
which we describe now. 



B. Bipartite purification strategy 

In the BEPP strategy, the parties A k , k = 1,2,. ..N 
wish to create a shared ensemble of iV-qubit graph states 
of high fidelity using a BEPP, where party A k holds the 
qubit corresponding to vèrtex a k - For each edge of the 
graph, one of the two parties connected by this edge pre- 
pares a connected graph-state pair, |G2,00) (equivalent 
to a Bell pair up to a local unitary) and sends one qubit 
of the pair to the other party through a noisy channel. 
(Alternatively, one could use a teleportation-based strat- 
egy: Alice distributes Bell pairs to the N — 1 other par- 
ties. The pairs are purified and then used to distribute 
the multipartite state that Alice has prepared locally.) 
The effect of the channels is given by Eq. I|13|) leading to 
states of fidelity F — q + and diagonal in the graph- 
state basis. The parties repeat the operation M times 
so that at the end M \E\ entangled pairs are distributed 
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between the different partners, where \E\ is the number 
of edges in the graph. The BEPP (reviewed in Section 
III D lli is then applied. This leads to a smaller ensemble 
of states given by a density matrix of the same form but 
with higher fidelity. Finally, the connection procedure 
described in Sec. III A6l is applied: Each party Ak merges 
together the |iV 0fc | qubits which will connect vèrtex a*, 
with its neighbors leading to the desired graph state. We 
call 



Y(F) 



final states # final states 



# initial states 



M 



the yield of the production of final states with fidelity F. 
To build up the desired multipartite state |G), we need 
one \G 2 ) pair for each edge of G. The number of edges 
for 1D cluster and GHZ states is (as for any tree graph) 
\E\ = N—l. Hence, the quantum communication cost is 
related to the yield by 



C 



F.G 



N-l 
Y(Fj- 



(27) 



The numerator is the number of channel uses (i. e. num- 
ber of transmitted qubits) required to distribute one 
state. This dependence on the size of the state prop- 
er ly reflects that for larger states, already the prepara- 
tion of the raw states is more costly. To allow for easier 
comparison with the yield, a figure that may fecl more 
familiar to the reader, we have plotted in all graphs the 



inverse communication cost C 
is proportional to the yield. 



F,G 



Y(F)/(N-1) which 



C. Multipartite purification strategy 

Alternatively, in the MEPP strategy, a central party, 
called Alice, creates M iV-qubit graph states locally. For 
each graph state, she keeps one qubit and sends the other 
N—l qubits through the channels to the N—l other par- 
ties. The resulting states are then purified using direct 
multipartite entanglement purification, i. e. the MEPP 
reviewed in Sec. III D 21 Hence, to distribute one state, 
we need N—l channel uses, the same as in the BEPP 
case. Thus, Eq. holds for MEPP, as well. 



D. Mixing of strategies 

Assume that the application of, say, m stcps of one of 
the protocols mentioned above reaches a final fidelity F\ 
with a communication cost of C\ , and application of m+ 1 
steps achieves fidelity Fi > F\ with communicatioin cost 
C 2 > Ci (i. e. Y2 < Yi). For a certain application, a fi- 
delity of F with Fi < F < F 2 is required, i. e. m steps are 
insufficient, but m+1 steps achieve a higher fidelity than 
desired at the cost of lower yield. In this case, one can 
find a compromise between the two strategies by mixing 
ensembles: 



Choosing an a € [0, 1], one prepares M raw states and 
then uses the first strategy on aM of them in order to 
gain aMYi states of fidelity Fi , and the second strategy 
on the remaining (1 — a)M states to obtain (1 — a)MY2 
states of fidelity F2 . Mixing these states gives an ensem- 
ble of fidelity 



aYíFi + (1 - a)Y 2 F 2 
aYx + {l- a)Y 2 



(28) 



with a a yield Y — ctY\ + (1 — a)Y 2 . This method allows 
one to obtain intermediate fidelities with a better yield. 
The communication cost mixes according to 



1 

C 



a 

cl 



1 



C 2 



(29) 



E. Intermediate strategies 

As a "compromise" between BEPP and MEPP, we 
shall also consider the following set of strategies: As- 
semble small states of N\ qubits, send them through 
the channels, purify them, and then use the connection 
scheme described in Sec. III A 61 to connect L of the N\- 
qubit states to one state with N = LN% — L + 1 qubits. 



IV. NUMERICAL SIMULATIONS 

A. Technique 

In genèric cases, an explicit numerical simulation of 
a quantum systems is intractable due to the exponen- 
tial growth of the Hilbert space with the number of in- 
volved partides or qubits. In our case, however, an effi- 
cient simulation is possible for two reasons: (i) All gates 
that are employed by the protocols are elements of the 
so-called Clifford group and hence, the Gottesman-Knill 
theorem applies, which allows for efficient simulations of 
pure state evolutions. (ii) The considered noise channels 
have Kraus representations that are diagonal in the Pauli 
basis. 

To explain (i), we start by reviewing the Gottesman- 
Knill theorem pill . It says that it is possible to sim- 
ulate so-called stabilizer circuits efficiently on a classi- 
cal computer. These are quantum circuits containing 
only preparation of computational basis states, opera- 
tions from the Clifford group, and measurements in the 
computational basis. The ./V-qubit Clifford group Cm is 
the group of those unitary operations that map Pauli 
operators onto Pauli operators under conjugation, i. e. 

C N := {Ue SU{2 N ) | UPUï eV N VP e V N } , 

V N :={±l,±i}·{l,(r x ,* y ,*,} 9N . (30) 

It happens to contain all the operations that we need for 
purifying, and hence, we can simulate the execution of 
the purification protocols described in Section IÏI Dl 
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Aaronson and Gottesman have given a fast algorithm 
which can perform such a simulation, and also supplied 
an implementation in the C programming language • 
We have used this software at the beginning of our stud- 
ies, but after realizing that its performance is not suffi- 
cient for our purposes, developed a new, faster algorithm, 
which is described elsewhere |24). 

The state represented in our simulator is always a pure 
state (Footnote |2q)- However, in entanglement purifi- 
cation, onc usually deals with mixed states, represented 
as density matrices. Nevertheless, due to the fulfillmcnt 
of condition (ii), we can get around this problem using a 
Monte Cario technique, which we describe now. 

To represent the ensembles of states we start with a 
high number Ni of qubits, typically several thousand 
times the number of qubits in the states to be purified. 
The qubits are initalized to a tensor product of \G, 0) 
states. Note that all these qubits can potentially get en- 
tangled, and hence have to be part of the same simulated 
quantum register. This would be prohibitive without a 
very efhcient algorithms for the stabilizer simulation. 

We then simulate all steps that are required to prepare 
Bell pairs or graph states, to purify them and to measure 
them. Depending on the measurement results, states are 
kept or discarded. Several iterations of the protocols are 
simulated. 

The transmission through the perfect channels amount 
to a simple relabeling: The program remembers the new 
site, where the qubit resides, as this indicates which 
qubits can be subject of joint operations. 

Simulating the channel noise is done by randomizing 
over many simulation runs as follows. The three noisy 
channels that we have considered, Eqs. (|1 lllHp . are sim- 
ulated using a pseudo-random number generator (RNG). 
Whcncver noise is to be applied onto a qubit, a random 
number between and 1 is generated, and if it is smaller 
than (1 — q) (the noise level), a x (a z ) is applied for bit- 
flip (phase-flip) noise. For depolarizing noise, the RNG 
is used again to obtain an integer between 1 and 4 which 
determines which of the operators 1, a x , a y , a z to apply. 

After the preparation of M initial states, m iterations 
of the protocol and for the BEPP case, connection of the 
purified pairs, Nf final states remain. The yield is then 
given by Y = This is, however, not a good estimate 
for the asymptotic yield in the limit of infinite ensembles 
for the following reason: If the number iVj_i of states at 
the beginning of purification step i is odd, we have to 
discard one state, because we can only deal with pairs of 
states. Hence, we better estimate the yield by 



Y 



N, 



N 



f 



(31) 



(with AI = No, Nf = N m , and |[ATJJ = N for even N, 
and [[ATJJ = N — 1 for odd N.) 

The fidelity can be determined by measuring the final 
states in the graph basis. This is because all the intended 
operations and the random noise operations map graph 
states onto graph states, so that all Nt final states are of 



the form \G,fi). The index fi can be determined as fol- 
lows: For each state, the graph state creation operation 
of Eq. ©, U{a.b}eE AZÍah) (which is Hermitian) is ap- 
plied again onto the state. If one then applies Hadamard 
gates on all qubits and measures in the a z basis, the 
measurement results spell out the index vector fi. As 
we intended to create \G, 0) states, we call the number of 
states for which we measured the number N g of "good" 
states and hence estimate the fidelity as 



N tot 

F = —f- ± 

A} ot 



/jVtot^tot _ jy-tot) 



(jVtot)3 



The superseript "tot" indicates that many runs of the 
simulations are made and that the numbers are the sums 
of the numbers in the individual runs. The uncertainty 
term follows from the expectation that, given a true fi- 
delity Ft, the number of good states F* ot output by 
the Monte Cario simulations after many runs is dis- 
tributed according to a binomial distribution with length 
Nj ot , hit probability Ft and hence Standard deviation 

Nj ot y/F T (l - F T ). Thus, is the estimate for F T with 
the given statistical uncertainty at la level. 

In the same way, the yield can be assigned an uncer- 
tainty, (Footnote [26j) 



Y 



N tot 

M tot 



± 



ÍNf t {M tot - Nj ot ) 

tot ^3 ' 



(M tot ) 



The la uncertainties are indicated by error bars in the 
plots. 



B. Extremal strategies 

We now present the results obtained for the two ex- 
tremal strategies described in Secs. IIII Bl and IÏII Cl with 
the following parameters: The distribution of the qubits 
is done through noisy channels and each step of proto- 
col requires imperfect two-qubit operations. The noise 
considered is depolarizing noise as defined in Eq. l(ï5j) 
with reliability p = 0.9 (10% noise) for the channels and 
pi = 0.99 (1% noise) for the local operations. We used the 
Monte-Carlo simulation method described in Sec. IIV Al to 
reach a precision on the fidelity varying from l%o to 1% 
depending on the size of the states and the number of 
iterations. 

For the MEPP case, one has to decide, which sequence 
of the sub-protocols PI and P2 to use. The alternat- 
ing sequence, P1-P2-P1-. . . , turns out to be not optimal 
in terms of yield and fidelity, neither for GHZ nor for 
cluster states. To find the optimum, one might hence 
consider to simulate, after each step, both sub-protocols, 
and then continue with the better one. Somewhat sur- 
prisingly, this leads to worse results (see Fig.|SJ). Thus, to 
find the optimal sequence of m protocol-steps, one would 
need to try all 2 m possibilities. As this is not practical, 
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FIG. 5: (Color online.) Comparison of the vàlues of inverse 
of communication cost and fidelity obtained after a number 
of steps varying from 1 to 5 for 3-qubit GHZ states. The red 
sòlid line stands for the BEPP strategy, the green dashed line 
for the alternating sequence of MEPP subprotocols beginning 
with PI and the blue small dashed line for the alternating 
sequence beginning with P2. 



we decided to stick with the alternating sequence, which 
turned out, though it is not optimal, to give very descent 
performance. For GHZ states, there is also a difference 
between the alternating sequences P1-P2-P1-. . . and P2- 
P1-P2-. . . , due to the asymmetry of the sets Va (con- 
taining only Alice's qubit) and Vb (containing the rest). 
Starting with PI works better, and this is what we use 
in all plots discussed in this section. 



1. GHZ states 

We start with the results obtained for GHZ states. We 
made our simulations for states of three to 10 qubits and 
a maximum number of steps varying from 5 to 7. As an 
example, Fig. [ü] shows the quantum communication cost 
as a function of desired fidelity for 5-qubit GHZ states. 
The data points are the outputs for 1 to 6 steps of the 
protocol. This plots allow us to determine, for a given 
fidelity, the strategy which will give the best yield (lowest 
communication cost). 

After 6 iterations, the increase in fidelity obtained by 
an additional step is smaller than the chosen precision 
of 1%. We therefore take this value as estimate of the 
maximum reachable fidelity. A comparison of the max- 
imal reachable fidelity for both strategies as a function 
of the number of parties (Fig. UJi) shows that the maxi- 
mal reachable fidelity is higher in the MEPP case for a 
number of parties strictly smaller than 10. In this case, 
there is a transition value of target fidelity from which on 
the MEPP strategy gives a better yield. We will refer to 
the value pair of fidelity and communication cost, where 
this transition happens, as the cross-over point. Fig. [S] 
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FIG. 6: (Color online.) Inverse of communication cost for 
different target fidelities for 5-qubit GHZ states and with p = 
.9 and pi = .99 (where p is the reliability defined in Eq. 1141 '). 
The data points are the outputs for 1, 2, . . . , 6 iterations of 
the protocol. The connecting lines are obtained by mixing 
ensembles of different fidelities according to Eq. The 
red sòlid line gives the obtained value in the MEPP and the 
green dashed line for the BEPP. The gain on fidelity from one 
step to the other becomes smaller at each step. From 6 to 7 
steps, the gain in fidelity is smaller than the uncertainty both 
in the BEPP and in the MEPP strategy. We consider this 
value as the maximal reachable fidelity. 



presents the yield as function of fidelity ior N = 3 and 
N = 10 as well as the cross-over points for intermediate 
number of parties. N — 9 is the highcst number of qubits 
for which there is a cross-over point. For higher number 
of parties, the BEPP strategy is always better. This is 
because of the fragilityjigainst noise of GHZ states for 
large particle numbers [271 l28j . The communication cost 
and fidelity of the cross-over points as function of the 
number of parties are presented in Fig. [ÜJl. 



2. Cluster states 

Next, we did simulations cluster states using the same 
parameter as for the GHZ states. The results are quite 
different. 

We made our simulations for states of three to fifteen 
qubits. In this range, as one can see in Fig. EU there is 
always a cross-over point. This is in stark contrast with 
the GHZ case. This main difference in behavior between 
this two kind of states is due to the much higher robust- 
ness of cluster states against noise [271 H3- Moreover, the 
range of target fidelity for which the multipartite strat- 
egy is the only one available increases with the number 
of parties as shown in Fig.[7p. In Fig.|ü^, we present the 
fidelity and communication cost of the cross-over point. 
Both vàlues decrease with the number of parties. This is 
due to the increasing cost of producing bigger and big- 
ger states and also to the fact that we consider here the 
global fidelity and not the LNE presented in Sec. III Cl 
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FIG. 7: (Color online.) Maximal reachable fidelity as function 
of N for (a) GHZ and (b) cluster states for the bipartite (green 
x ) and multipartite (red bars) strategies with reliabilities (cf. 
Eq. CU) P = 0-9 f° r channel transmission and pi — 0.99 for 
local operations. The final fidelity is estimated as follows: 
For a given number of parties, we iterated the protocol as 
long as we obtained an increase of fidelity larger than the 
uncertainty (typically 1%). We took the last value as maximal 
fidelity and assigned its uncertainty to the maximal reachable 
fidelity. The green crosses give the vàlues in the bipartite 
case while the red bars give it for the multipartite case. One 
sees here the main difference in behavior between GHZ and 
cluster states. In the first case, there is a range where the 
multipartite strategy is better than the bipartite one for a 
number of parties strictly smaller than 10. For more parties, 
the multipartite protocol fails because of the fragility of GHZ 
states against noise. On the other hand, the robustness of 
cluster states allow us to purify them even for a large number 
of parties. The range of fidelity where MEPP is superior 
increases with the number of parties. 
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FIG. 8: (Color online.) Inverse of communication cost for 
different target fidelities of GHZ states of 3 to 10 qubits with 
alteration probabilities (as in Fig. 1 — p = 0.1 and 1 — pi = 
0.01. The dashed green line stands for 3-qubit GHZ states and 
BEPP strategy, the red sòlid line for 3-qubit GHZ states and 
MEPP strategy, the blue small-dashed line for 10-qubit GHZ 
states and BEPP strategy, the pink dotted line for 10-qubit 
GHZ states and MEPP strategy. The blue squares give the 
cross-over points, i.e. the fidelity where MEPP becomes more 
efficient than BEPP, for N = 3,5, 6, 7, 8 and 9. For N = 3 
and N = 10, the purification curves are plotted as well. For 
N = 3, they cross at the corresponding blue squarc. 
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FIG. 9: (Color online.) Inverse of communication cost (green 
x) and fidelity (red +) of the cross-over depending on the 
number of parties for (a) GHZ and (b) cluster states. This 
vàlues are obtained by using the Monte Cario method and 
are therefore submitted to errors. The cross-over indicates 
the range of target fidelity from which up the MEPP strategy 
is more efficient than the BEPP strategy. Note the log scale 
for the inverse of communication cost. In the GHZ case, there 
is no cross-over for more than 9 parties. 
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FIG. 10: (Color online.) Inverse of communication cost 
for different target fidelities for 3 (red sòlid line for MEPP 
and green dashed line for BEPP) and 15 (pink dotted line 
for MEPP and blue small dashed line for BEPP) qubit 
cluster states. The data points are the outputs for 1, 2, 
3, . . . iterations of the protocol. The intermediate points are 
obtained by mixing ensembles of different fidelities. For more 
than 6 steps, the difference between the reached fidelity and 
the maximum reachable fidelity is smaller than the uncer- 
tainty. For any number of parties, the curves representing 
the two strategies cross over. The disks give this cross-over 
for N = 3, 4, 5, 6, 7, 8, 9, 10, 15. (That one curve seems to "go 
back" is just an artifact of the statistical inaccuracies of the 
Monte Cario method.) 
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Un prepare an n-qubit cluster or GHZ state [JJJ as well as Fíg.HTl show a regauging of the fidelity axis 

B2 prepare a Bell pair to the SO-called local noise equivalent (LNE) described in 

S send the states through the channels g gc W\J^\ 

PI apply multipartite purification protocol PI 

P2 apply multipartite purification protocol P2 

Pb apply bipartite protocol 

Cl connect í states to a larger one 

TABLE I: Legend for the instruction strings in Figs. 1111 and 

m 



C. Intermediate strategies 



Since switching from BEPP to MEPP can result in 
such striking differences in yield, one might expect that, 
especially near the break-even point, certain intermediate 
strategies, mixing characteristics of BEPP and MEPP, 
might perform even better. After all, in the BEPP case, 
one purifies small states (with only 2 qubits) and then 
connects them, while in the MEPP scenario, the states 
are first connected to large units, which are then purified. 
One can also connect pairs to states of intermediate size, 
purify these, connect them to the dcsired full size, and 
perhaps purify again. This can be seen e. g. in Fig. llll In 
this figure, we have simulated many different strategies 
which are described in short by instruction strings which 
are processed from left to right and teli the software in 
which order which preparations, transmissions, connec- 
tions or purihcations should be simulated (cf. Table[IJ|. 

It can be seen that for low fidelities and high yields 
(left side of the plots), the BEPP case is best, as already 
seen above, and for high fidelities and low yields (right 
margin of the plots), MEPP catches up. In the middle 
rcgion, one may indeed increase the performance by hrst 
preparing small states of, say, 4, 5 or 7 qubits, purifying 
them, and then connecting them to the desired 13-qubit 
state. (Do not get confused by the appearance of "M13- 
S" at the left margin. This looks like MEPP, but is not, 
as it contains no purification at all. Also note that there 
is a subtle difference between using the BEPP protocol 
(denoted "B2-S-Pb-. ..") and using the MEPP protocol 
on the \G 2 ) state (denoted "M2-S-P1-. . ." or "M2-S-P2- 
..."), with the former performing better.) 

Of course, only discrete ways of assembling the desired 
states from equal smaller states are available. Recali that 
connecting L states of n qubits will give a state of 

N = Ln-(L-í) (32) 

qubits because (L — 1) qubits have to be measured in the 
connection process. In the plots, we have taken all pos- 
sible vàlues of L for the given state size n and calculated 
data points for the corresponding strategies with up to 
four purification steps. The blue curve in the plot marks 
the optimum that can be achieved using theses strategies, 
and mixing them as described in Section FlII Dl 

To demonstrate the efficiency of our procedure, we 
also considered purification of cluster states of 31 qubits 
(Fig. 1121) . In order to allow for easier comparison, Fig. 
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FIG. 11: (Color online). Examples for the use of intermediate purification strategies, here for 13-qubit (a) cluster and (b) GHZ 
states. Plotted is the inverse of communication cost as function of final fidelity. See Table[I]for the meaning of the instruction 
strings. The blue curve marks the maximal achievable yield for a given desired fidelity and is obtained by connecting the 
optimal strategies with curves according to Eq. 11281 . Noise levels are (1 — q) = 0.1 for the channels and (1 — qi) = 0.01 for 
local operations. In (a), one can -following the dark blue line- see well, how for small target fidelity BEPP ("M2-S-. . .") gives 
the best yield, while for high fidelities (F > 0.9), distributing larger and larger states becomes advantageous. (For even higher 
fidelities, one expects the full MEPP strategy, i. e. "M13-S-. . .", to appear on the blue curve. However, this will happen at 
communication costs larger than the scales shown on the plot, which ends with "M13-S-P1-P2-P1" , i. e. MEPP with only three 
purification steps.) In (b), the picture is not as clear, at the GHZ states already start to deteriorate under the given level of 
local noise. 
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FIG. 12: (Color online). Production of 31-qubit cluster states, using intermediate strategies. The "instruction strings" are 
explained in Table Data points with the same number of purification steps are plotted in the same color. Note how the 
distribution of initially larger states becomes advantegeous for higher target fidelities. Noise levels are (1 — q) = 0.1 for the 
channels and (1 — qi) = 0.01 for local operations. (Note also that the the data points at the low end of the plot have few 
purification steps: Those steps beginning with "M16" or "M31" that should appear on the curve of aptimal strategies are again, 
as in Fig. 1111 beyond the range of the plot.) 



V. ANALYTICAL TREATMENT FOR A 
SIMPLIFIED MODEL 

For a better understanding of the numerical results, 
we now develop an analytical treatment for both BEPP 
and MEPP. To makc this task feasible we have to restrict 
ourselves to a simplified noise model. We only consider 
GHZ states. 

As before, we define two sets Va and Vb corresponding 
to the bi-coloration of the graph. Va is the set contain- 
ing only one qubit, namely the central vèrtex which is 
connected to all the others, and Vb contains the rest. 

In the toy models presented below, the central party, 



called Alice wants to share an iV-qubit GHZ state with 
(JV — 1) partners. Depending on the strategy, the initial 
states are either Bell pairs or GHZ states, which are noisy 
due to the transmission through the channels. First, in 
Subsection IV Al we consider local operations to be per- 
fcct. We will see that this fails to reproduce features seen 
in the numerical results. Hence, we extend our model, in 
Section lVBI such that it incorporates local noise. 
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A. Perfect local operations 

To start, we assume to local operations to be per- 
fect. Of the channels we considered in Eqs. I|11I13|I . only 
bit-flip channels and phase-flip channels allow for a sim- 
ple analytical treatment. We present the calculation for 
phase-flip channels. The calculation and the results for 
bit-flip channels are very similar. We have hence not 
included them in the paper. 



1. BEPP strategy 

Following the BEPP scenario described in Sec. Ull Bl 
Alice sends one qubit of each entangled pair p = 
|G 2 ; 0, 0) (G 2 ; 0, 0|of her initial ensemble through the 
channel to party B^, obtaining 



k = 1, . . . , N — 1, in order to obtain the iV-qubit GHZ 
state. 

Note that the qubits that Alice connects have not been 
sent through channels and are hence unaffected by chan- 
nel noise. Thus, it does not matter whether we first apply 
the superoperator for the channel noise and then the one 
for the local noise due to the connection process, or vice 
versa. This means that the final state is obtained by ap- 
plying noise on all qubits of the GHZ state that do not 
belong to Alice. This leads to a fidelity 



F 



BEPP 



(iV,ç,m) = / BEPP (ç,m) 



and (as the channels are used N — 1 times to create one 
-/V-qubit GHZ state) a quantum communication cost 



C 



BEPP 



N - 1 



Y BEPP (q,m)' 



p = q\G 2 ;Q, 0)(G 2 ; 0, 0| + (1 - ç)|G 2 ;0, 1)(G 2 ; 0, 1|. 

She then applies the BEPP. The state of the pairs that 
are kept after one step is given by (see Eq. ITÜ|l 



|G 2 ;0,0)(G 2 ;0,0| 



q 2 + {l-q) : 



+ 9 (1 9)2 |G 2 ;0,1)(G 2 ;0,1|. (33) 
T + (1 - 9) 

In this step, the probability of keeping the source state 
after the measurement of the other state is given by 
£; BEPP (<7, 1) = q 2 + (1 — q) 2 (probability of having same 
measurement outcomes). We denote by / BEPP (g,m) the 
fidelity after m steps. The quantity k(q, m) is called 
success probability in step va. Note that the ratios of 
the ensemble sizes after and before the step is given by 
k(q, m)/2 as one half of the states are measured and dis- 
carded. One obtains the total yield y BEPP after m steps 
by multiplying these ratios for the individual steps. By 
iterating the protocol over m steps one finds: 



/ BEPP ( ?J m) 
fc BEPP ( g ,m) 

y BEPP ( g ,™) 



q^ + (1 - qf m ' 
q 2m + (1 - qf m 



+ (l-9) 2 



(34) 
(35) 



n 

i=l 



k(q,i) 



q 2m + (1 - qf 



2 m n™i (q^ + il-q) 2 



(36) 



After the bipartite purification, Alice connects (N — 1) 
pairs to produce an TV-qubit GHZ state. To connect two 
pairs, she applies a controlled phase gate (Eq. ©) fol- 
lowed by a a y measurement on one of the two qubits 
just connected (cf. Fig. |3J). This procedure is repeated 
(N — 1) times between different pairs of parties (A, Bk), 



2. MEPP strategy 

In the MEPP setting, Alice prepares an 7V-qubit GHZ 
state locally and distributes it through depolarizing chan- 
nels to her (N — 1) partners. We then have the state 



\a=2 



(a) 



|G*;0,0)(G*;0,0| 



(37) 



where is the phase-flip channel dcfined in Eq. Ullfl . 
We shall from now on suppress the symbol G* which 
indicates the iV-vertices star graph of Fig. [21 

We shall see that all states that we encounter have the 
form 



r, 



(m) 



I) 



0)<0| 



N 
ii=2 



1 

.010. 



■0)(0,. 



-r. 



N 

íl ,12=2 

h<Í2 



where denotes the coefficient in front of the terms 
with j entries "1" after the mth step of the purification 
protocol. These states are diagonal in the graph state 
basis and symmetric w. r. t. permutations of the qubits 
in set Vb ■ They are hence characterized by only N coef- 



|0,. 



1 

.010. 



1 

.010. 



■><o, 



1 

.010. 



1 

.010. 



.0|- 



1 

.010. 



tJSxI0.11. 



1)<0,11...1|, (38) 



„( m ) 



„( m ) 



ficients rg ,...,/ jv-x- 

We start by carrying out the application of the super- 
operator in Eq. H37[). Indeed, one obtains a mixture of 
the form l|38|l with coeficients 



,(°» 



<1 



N-l-j 



(i - q y 



(39) 



As only set Vb is affected by the noise, subprotocol 
P2 is sufficient to purify the state. Following 0, one 
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sees that after each step of the subprotocol, the state is 
changed such that each coeficient becomes proportional 
to the square of its former value, i. e. 



-("») = 

3 



„(«»-!) 



l~,i=0 \ i ) 



(m-1) 



2 ' 



(40) 



Inserting Eq. (JSÜJ, one gets for the first step (using the 
binomial theorem) 

(l) _ q 2(N-l- j){1 _ q) 2 j 

" j [f + (1 - q)*]»- 1 ' 
and iterating the formula, one finds 

[q** + (1 - qf m f- 1 ' 
The fidelity of the state at step m can now be read off: 



F 



MEPP 



(TV, q, m) = r [ Q 



(m) 



q 2m + (1 - q) 2 



N-l 



Note that this is the same expression as we got before for 
the BEPP: F MEPP (N, q, m) = F BEPP (N, q, m). 

To calculate the yield, we need the success probabil- 
ity k MEPP (N, q,m) that a state is kept. Using a similar 
argument as we did for Eq. <|4U|I we finds 



N-l 



fc MEPP (/V,ç,m)=]T 



N — 1 



(m-l) 



q 2m + (1 - qf 



N-l 



V [q 2 - 1 +(l-q) 
From this, we can find the yield as before 

Y MEPP (N, q,m) = H - 

_ n - m ( q 2m + (i - q? m 



-n 2 



N-l 



vn^T 1 [q^ + a-qr]) 

Comparing to the BEPP case, one sees that 



(41) 



Y 



MEPP 



(N,q, m) 



_ nm(Af-2) [yBEPP 



[Y^ PP (N,q,m)] 



N-l 



(42) 



3. Conclusion 



In the particular case of dephasing channels and per- 
fect local operations, both strategies lead to the same 
fidelity after iterating the protocol the same number of 
steps. However, they differ in the communication cost. 



As one sees from Eq. (H2*|l , the yield of the MEPP strat- 
egy is always lower (and the communication cost hence 
larger) . This fact can be explained from the higher prob- 
ability of throwing away states at each step, which even 
increases further with the number of parties. 

We have also done analytical calculations for bit-flip 
channels (Eq. l|l·lf) ) and numerical simulations for depo- 
larizing channels (Eq. Ijl^f) ) (always in case of perfect 
local operations) and found a similar behavior. 

In order to see regions where MEPP is superior as we 
did with the Monte Cario simulations, it is hence neces- 
sary to give up the simplification of assuming noiseless 
local operations. 



B. Imperfect local operations 

If local operations are not assumed to be perfect, re- 
sults are quite different. We again consider GHZ states 
of arbitrary size, to be purified with the BEPP or MEPP 
strategy. 

For the noise, we define a model which is simple enough 
to allow for analytical calculations but still shows the gen- 
eral features obtained numerically (especially it shows 
cross-over points): The channel through which Alice 
sends qubit to the other parties is the phase-flip channel 
of Eq. Ijll|l with alteration probability (1 — q). The im- 
perfection of the local gates are modeled by bit-flip noise 
(Eq. jXHO , f° r Alice's operations, and phase-flip noise 
for all other parties, always with alteration probability 

(!-«)• 



1. MEPP strategy 

As before, Alice prepares perfect iV-qubit GHZ state 
and distributes them using the channels. We get the 
same initial state as before, again of the form l|38(l 
with coefficients as in Eq. (I39f) . We shall see that, again, 
the form (|38|l is preserved by the purification steps even 
though they are now assumed to be noisy. 

The vàlues of the are changed according to a linear 
map: 

We shall construct this map in two steps. First, we 
see, how the phase-flip noise of the local gates acting on 
the qubits in Vb, 



N 



(43) 



changes the coefficients and denote the map correspond- 
ing to this action by A: 



N-l 



k=0 
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Then, we consider the action of the bit-flip noise on Al- 
ice's qubit get the full map A. 

For the first step, we call a state |G* , fi) a fc-state, if \x 
starts with a (for the central qubit in Va) and contains, 
within the indices corresponding to Vb, k entries "1" 
and (N — 1 — k) entries "0" . We can now calculate the 
probability Pj^k that the superoperator 1|43[) changes a 
pure fc-state to any j-state: Say, s of the k entries "1" 
are flipped to "0" . Then s = j — k + soï the (N—l — k) 
entries "0" have to be flipped to "1". Hence, 



£ I (i-*)-«í 



< ^"^^(l-flïW- 1 -*- 1 . ( 44 ) 



There are ( . 1 ) j-states and ( k 1 ) fc-states, and so 



\ k 

(N-l\ 

v i 



~ 'tN-í\Pi^ k i 



which can also be written in terms of Gaufi's hypergeo- 
metric function F as 



Xjk ~ fc!(j + fc)f 



\ F I ./' ' ; -k, i - N + 1, i + k + lj x 

x(l- qi y- k q?- X -i +k (45) 

Now, we can do the second step and apply the noise 
for the imperfection of Alice's local gates, modeled by 
£x\qi). We get, due to Eq. ©, 



N-l 



qi r' 3 + (1 - qi )r^_, = J2 A. 



krk, 



k=0 



where 



N-i\ k 



A 



V j ) s=0 v 

N-l-k 
j -k + s 
N-l-k 
3-s . 



N-l-j+k-2s+l^ _ qi y-k+2s + 



1l 



N-l-j-k+2s 



il _ qi)j+ k-2s+l 



(46) 



The fidelity and the yield corresponding to one step of 
protocol can then be calculated: 



f(N,q,qi,m) = 



2^k=o { k ) 



„("») 



E,=o A 0J r 



(m-l) 



EiV-1 /JV-1\ » 
fc=0 V fc J ^2^j=o Aí 



Jm-l) 
3 



As before, the denominator of the previous expres- 
sion is the success probability k(N,q,qi,m) of the 
step to (Footnote 29]), and we get the total yield 
Y MEFP (N,q,qi,m) by multiplying up the factors k/2 of 
all m steps: 



Y^(N,q, qi ,m) = l[ 



k{N,q, q u i) 



(For the corresponding quantum communication cost, 
Eq. ((23) has to be used.) 



BEPP strategy 



Next, we find an analytical treatment for the BEPP 
strategy. In order to facilitate the calculation, we will 
consider the BEPP as a special case of the MEPP. We 
first show why this is possible without changing the re- 
sults: 

Considering the restricted noise model presented in 
this Section, the state can always be written as a contri- 
bution of |G 2 ;0,0)(G 2 ;0,0| and |G 2 ; 0, 1) (G 2 ; 0, 1| only. 
It can hence be purified using only subprotocol P2. In 
addition, the only difference between the BEPP protocol 
as described in Sec. 111 DU and subprotocol P2 is the ex- 
change of the states |G 2 ; 1, 0) and |G 2 ; 1, 1) in the former, 
and not in the latter. As these states have no contribu- 
tion in the present model, the two protocols give identical 
results. 

The results obtained in the last section can be 
used to calculate the fidelity f(q,qi,m) and the yield 
Y (q, qi,m) before connection. After purification, 
Alice connects qubits a%, . . . ,a/y-i of (N — 1) pairs 
described by the states pk = ç|G 2 )/ £ (G2| + (1 — 



q)ví bk) \G 2 )k(G 2 \aí bk} 
pairs is given by 



respectively. The joint state of the 



JV-l 



JV-l 



/' = '/' ' I ^9 |G 2 ) fc (G 2 | 

N-l /JV-l \ 

* } ' «§\G 2 ) k (G 2 \ 4 6i) + ··- 
í=i \fe=i / 

/N—X \ /N-l \ /N-l 



-(i- «)í"- a 5>;w 



The connection is performed using the procedure de- 
scribed in Sec. III A 61 As the noise contains only a z op- 
erators, it commutes with the projectors and correction 
operators that are used in this procedure. The state after 
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projection is given by 

P 2 p = q N - 1 \G*)(G*\ 



N 



b=2 

Hi-if- 1 fn^ 6) ) i^x^i (fi^A- 

It follows that after the connection, the fidelity is given 
by 

F BEPP = f(q,qi,m) N - 1 
and the quantum communication cost by 



C 



BEPP 



N- 1 



Y MEPP (2,q, qi ,m) 



3. Discussion 

The results obtained the way just explained are pre- 
sented in Figs. El d and [13 always for MEPP and 
BEPP. 

The cross-over points, i. e. the fidelity vàlues (and 
corresponding costs) above which MEPP performs bet- 
ter than BEPP, are plotted as blue squarees in Fig. El 
for up to 70 parties. We have also plotted the fidelity- 
cost function for selected numbers TV of parties, which 
cross in their respective blue squares. The other cross 
points were determined in the same way by observing 
where the MEPP and the BEPP curve intersect. Com- 
pare this figure with Fig. |S] Our analytical toy model 
could reproduce the appearence of cross points which we 
had already discussed in Section IÏVB II an essential fea- 
ture observed for the more general noise model. It does 
not, however, reproduce the fact that a cross-over ceases 
to appear above a certain number (here: 9) of parties. 
This fact is due to the particular kind of noise of our 
toy model under which GHZ states appear less fragile 
than under depolarizing noise so that the break-down of 
MEPP for large states does not happen. 

We can use our analytic model to explore the param- 
eter space more thoroughly. For instance, one might be 
interested how the positions of the cross points changc 
if the local noise is increased. This is shown in Fig. ITT1 
where the right-most curve is the same as the disks in 
Fig. El and the others are for higher local noise levels. 
Observe how the effect of local noise depends more and 
more on the state size as its level approaches the order 
of magnitude of the channel noise. 

The vertical tails of the curves in Fig. ^] already allow 
to easily read off the maximum reachable fidelity, which 
is plotted in Fig. El There, the advantage of MEPP 
over BEPP increases with the number of parties. This 
effect can also be seen in the numerical calculations for 
depolarizing noise (Fig. UJi). In the latter case, it is, 
however, soon overcome by the competing effect of the 
break-down of MEPP under realistic (depolarizing) noise. 
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FIG. 13: (Color online.) Inverse of communication cost as 
function of final fidelity for the simplified noise model used in 
section IV Bl Analytical calculation for GHZ states of differ- 
ent number of qubits N varying from 5 to 70, with alteration 
probability for the channel and local noise of (1 — q) = 0.1 
and (1 — qi) = 0.05 respectively. The green dashed lines stand 
for MEPP strategy while the red sòlid lines stand for BEPP 
strategy. The blue circles give the crossing points for all num- 
ber of parties between 5 and 70. 
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FIG 14: (Color online.) Analytical results for different num- 
ber of parties and different amount of local noise. Each curve 
gives the yield as function of fidelity for the cross-over for a 
given alteration probability (1 — qi) (See Eq. lli U and jH} ). 
This parameter varies from qi = 0.93 (left curve) to qi = 0.99 
(right curve). 



C. Testing the numèrics 

The analytical formulas are also very useful for verify- 
ing the code of our numerical calculations. Switching the 
programs from depolarizing noise to the simplified noise 
considered here is a trivial alteration. We find that the 
numerical results agree well with the analytics, see Fig. 
El t° El This fact makes us confident in the correetness 
of our codes. 
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FIG. 15: (Color online.) Maximal reachable fidelity F max 
plotted against the number of parties for the simplified model 
described in Sec. IV Bl applied to GHZ states. The alteration 
probability for the channels and the local noise are given by 
(1 — q) — 0.1 and (1 — qi) = 0.01 respectively. The results 
were obtained analytically. 
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FIG. 16: (Color online.) Testing the numèrics: We switched 
the programs that were used to calculate the results of section 
II V Bl íprogram C) and section IIV Cl fprogram S) to the simpli- 
fied noise model of section lVBI The plot shows the inverse of 
communication cost as function of the final fidelity. The red 
symbols (+) stand for the analytical results while the green 
(x) and the blue (*) symbols stand for the output of pro- 
gram C and program S respectively. The error bars stand for 
la errors. A comparison with the derived analytical formu- 
las shows satisfactory agreement. The calculation was done 
for MEPP of GHZ-states with 10 qubits (b) at noise levels 
q = 0.9 and q t = 0.95. 



VI. SUMMARY AND CONCLUSIONS 

In this article, we have investigated the quantum com- 
munication cost of preparing a class of multipartite en- 
tangled states with high fidelity. The presence of noisy 
quantum channels and imperfect local control operations 
rcquires the usage of error correction or -in our case- 
entanglement purification schemes to achieve this aim. 



Ü 



o 
o 

E 
E 
o 
o 



0.1 



0.01 



I I 


BPP, numerical i 1 • 




BPP, analytical 






; 


M 




H< " 




*• 

1 1 1 1 



0.93 0.94 0.95 



0.96 
F 



0.97 0.98 0.99 



fin 



FIG. 17: (Color online.) Testing the numèrics obtained for the 
BEPP strategy. Noisy entangled pairs, arising from sending 
one of the qubits through a depolarizing channel, are purified 
using the BEPP. The plot shows the inverse of communication 
cost as function of final fidelity. The black crosses give the 
exact vàlues while the red bars give the numerical results of 
the Monte-Carlo simulation. 



We have considered various strategies to generate these 
high-fidelity states and have established in this way up- 
per bounds on the quantum communication cost. The 
optimal strategy strongly depends on the error param- 
eters for channels and local control operations, and on 
the desired target fidelity. For a simple error model and 
the generation of GHZ states based on various strategies, 
we have obtained analytic results that allow us to com- 
pare these strategies. Numerical simulations for genèric 
error models, based on Monte Cario simulation, show es- 
sentially the same features as observed in the simplified 
model. The simulation makes use of a recently developed 
method that allows one to efficiently simulate the evolu- 
tion of stabilizer states (or graph states) under Clifford 
operations on a classical computer |23l |24| . We have also 
applied this method to investigate not only the genera- 
tion of GHZ states but also of other types of multipartite 
entangled states, e.g. cluster states. 

We find that for high target fideüties, strategies based 
on multipartite entanglement purification generally per- 
forin better than strategies based on bipartite purifica- 
tion. For low target fidelities, strategies based on bi- 
partite purification have a higher efl·lciency, leading to 
smallcr communication cost. 

We believe that the generation of high-fidelity multi- 
partite entangled states is of significant importance in the 
context of (distributed) quantum information processing. 
Such multipartite entangled states represent resources, 
e. g. for measurement-based quantum computation, con- 
ference key agreement and secret sharing schemes, and 
may be used for other security tasks. Our investigation 
takes both channel noise and noisy apparatus into ac- 
count. We could show that the choice of a proper strategy 
not only allows one to significantly reduce the quantum 
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communication cost, but to reach fidelity of target state 
that arc not accessible otherwise. 
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